NONLINEAR ACCELERATOR PROBLEMS VIA WAVELETS: 
6. REPRESENTATIONS AND QUASICLASSICS VIA FWT 
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In this series of eight papers we present the applications of 
methods from wavelet analysis to polynomial approxima- 
tions for a number of accelerator physics problems. In this 
part we consider application of FWT to metaplectic repre- 
sentation(quantum and chaotical problems) and quasiclas- 
sics. 

1 INTRODUCTION 

This is the sixth part of our eight presentations in which 
we consider applications of methods from wavelet anal- 
ysis to nonlinear accelerator physics problems. This is a 
continuation of our results from [l]-[8], in which we con- 
sidered the applications of a number of analytical methods 
from nonlinear (local) Fourier analysis, or wavelet analy- 
sis, to nonlinear accelerator physics problems both general 
and with additional structures (Hamiltonian, symplectic or 
quasicomplex), chaotic, quasiclassical, quantum. Wavelet 
analysis is a relatively novel set of mathematical methods, 
which gives us a possibility to work with well-localized 
bases in functional spaces and with the general type of 
operators (differential, integral, pseudodifferential) in such 
bases. In contrast with parts 1-4 in parts 5-8 we try to 
take into account before using power analytical approaches 
underlying algebraical, geometrical, topological structures 
related to kinematical, dynamical and hidden symmetry of 
physical problems. In part 2 according to the orbit method 
and by using construction from the geometric quantization 
theory we construct the symplectic and Poisson structures 
associated with generalized wavelets by using metaplectic 
structure. In part 3 we consider applications of very useful 
fast wavelet transform technique (FWT) (part 4) to calcu- 
lations in quasiclassical evolution dynamics. This method 
gives maximally sparse representation of (differential) op- 
erator that allows us to take into account contribution from 
each level of resolution. 

2 METAPLECTIC GROUP AND 
REPRESENTATIONS 

Let Sp(n) be symplectic group, Mp(n) be its unique two- 
fold covering - metaplectic group [9]. Let V be a sym- 
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plectic vector space with symplectic form ( , ), then R © V 
is nilpotent Lie algebra - Heisenberg algebra: [R, V] — 

0, [v,w] = (v,w) S R, [V,V] = R. Sp{V) is a group 
of automorphisms of Heisenberg algebra. Let N be a group 
with Lie algebra R ffi V, i.e. Heisenberg group. By Stone- 
von Neumann theorem Heisenberg group has unique irre- 
ducible unitary representation in which 1 i— ► i. Let us also 
consider the projective representation of symplectic group 
Sp(V): U gi U g2 = c(g 1: g 2 ) ■ U gig2 , where c is a map: 
Sp(V) x Sp(V) — > S 1 , i.e. c is S' 1 -cocycle. But this rep- 
resentation is unitary representation of universal covering, 

1. e. metaplectic group Mp(V). We give this representa- 
tion without Stone-von Neumann theorem. Consider a new 
group F = N 1 XI Mp(V), txi is semidirect product (we 
consider instead of N = R © V the N' = S 1 x V, S 1 = 
(R/2wZ)). Let V* be dual to V, G(V*) be automorphism 
group of y*.Then F is subgroup of G(V*), which consists 
of elements, which acts on V* by affine transformations. 
This is the key point! Let qi, ...,q n ;pi, ...,p n be symplec- 
tic basis in V, a — pdq — J^Pidqi and da be symplectic 
form on V*. Let M be fixed affine polarization, then for 
a £ F the map a i— > 6 a gives unitary representation of G: 
9 a : H(M) — > H(M). Explicitly we have for representa- 
tion of N on H(M): (@ g f)*(x) = e-* qx f(x), Q p f(x) = 
f(x — p) The representation of N on H(M) is irreducible. 
Let A q , A p be infinitesimal operators of this representation 



A n 



lim - [6 



[Q- tg -l], 



then Agf(x) = i(qx)f(x), A p f(x) = 



^dx~ {x) 



Now we give the representation of infinitesimal basic ele- 
ments. Lie algebra of the group F is the algebra of all (non- 
homogeneous) quadratic polynomials of (p,q) relatively 
Poisson bracket (PB). The basis of this algebra consists of 
elements l,qi,...,q n , pi,...,p n , qiqj,qiPj, PiPj, i,j = 
l,...,n, i<j, 



df dg 
^ dpj dqi 



dldg_ 

dq l dpi 



PBis {f,g} 

and {1,3} = for all g, 
{Pi ,Qj} = % , {Pilj , Ik } = S ik q 3 , 
{PiQj,Pk} = SjkPi, {PiPj,Pk} = 0, 
{PiPj , Qk } = SikP 3 + 5jkPi , {qiqj , qk } = 0, 

{q%qj,Pk} = -SikQj - 5 jk qi 



so, we have the representation of basic elements / i— ► Af 

1 i-> i, q k i-> iarfc, 



<9 



i a* 



i dx k dx l 



,QkQl 



ix k x l 



This gives the structure of the Poisson manifolds to rep- 
resentation of any (nilpotent) algebra or in other words to 
continuous wavelet transform. According to this approach 
we can construct by using methods of geometric quantiza- 
tion theory many "symplectic wavelet constructions" with 
corresponding symplectic or Poisson structure on it. Then 
we may produce symplectic invariant wavelet calculations 
for PB or commutators which we may use in quantization 
procedure or in chaotic dynamics (part 8) via operator rep- 
resentation from section 4. 

3 QUASICLASSICAL EVOLUTION 



where a = (pt\, . . . , a n ) is a multi-index, \a\ = cti + . . . + 
oi n , D x — —iUd x . So, evolution (1) for symbol bt(x, £; %) 
is 



h = {H,b t } + ^j E 

\a\+P\=j>l 

W{d?HD%).{dh t D<* x H). 



(2) 



At h = this equation transforms to classical Liouville 
equation 

b t = {H,b t }. (3) 

Equation (2) plays a key role in many quantum (semiclas- 
sical) problem. We note only the problem of relation be- 
tween quantum and classical evolutions or how long the 
evolution of the quantum observables is determined by the 
corresponding classical one [9]. Our approach to solution 
of systems (2), (3) is based on our technique from [l]-[8] 
and very useful linear parametrization for differential oper- 
ators which we present in the next section. 



Let us consider classical and quantum dynamics in phase 
space ft = R 2m with coordinates (x, £) and generated by 
Hamiltonian H(x,£) G C°°(ft;i?). If $^ : ft — > ft is 
(classical) flow then time evolution of any bounded clas- 
sical observable or symbol b(x,£) G C°°(ft, R) is given 
by b t (x,£) = b(<f>?(x,£)). Let H = Op w (H) and 
B = Op w (b) are the self-adjoint operators or quantum 
observables in L 2 (R n ), representing the Weyl quantization 
of the symbols H, b [9] 



1 



(Bu)(x) = 

(27r7l) r - j R 2„ 
e i < x - yU> / n u{y)dyd^, 



y 



where u G S(R n ) and B t = e iHt/h Be -iHt/h be ^ 
Heisenberg observable or quantum evolution of the observ- 
able B under unitary group generated by H. B t solves 
the Heisenberg equation of motion B t = {i/h) [H, B t ] . Let 
bt(x,£;U) is a symbol of B t then we have the following 
equation for it 

b t = {HA}M, (1) 

with the initial condition bo(x,£, Ti) = b(x,£). Here 
{/, g}M{x,Q is the Moyal brackets of the observables 
f,g€ C°°(R 2n ), {f,g} M (x,0 = ftg - «?«/, where f$g 
is the symbol of the operator product and is presented by 
the composition of the symbols /, g 



(/tt<7)(*,0 



1 



(2tt?j.)«/ 2 



i<r.p> /h+i<u,T>/h 



■f(x + lj, p + £) ■ g(x + r,r + £,)dpdrdrduj . 

For our problems it is useful that {/, <?}m admits the formal 
expansion in powers of fi: 

{/, sM*,0~ {/.»} + 2- J '- 

E (- 1 ) 1 ' 3 ' • We ' i^gD^f), 

\a+0\=j>l 



4 FAST WAVELET TRANSFORM FOR 
DIFFERENTIAL OPERATORS 

Let us consider multiresolution representation . . . C V<i C 
V\ C Vo C V-i C V-2 ■ ■ ■ (see our other papers from 
this series for details of wavelet machinery). Let T be an 
operator T : L 2 (R) -> L 2 (R), with the kernel K(x,y) 
and Pj : L 2 (R) — > Vj (j G Z) is projection opera- 
tors on the subspace Vj corresponding to j level of res- 



olution: (Pjf)(x 



Qj = Pj-i — Pj is the projection operator on the sub- 
space Wj then we have the following "microscopic or tele- 
scopic" representation of operator T which takes into ac- 
count contributions from each level of resolution from dif- 
ferent scales starting with coarsest and ending to finest 
scales [10]: T = V , X U) ,IU , + Q 3 TP 3 + PjTQj). We 
remember that this is a result of presence of affine group in- 
side this construction. The non-standard form of operator 
representation [10] is a representation of an operator T as 
a chain of triples T = { Aj , Bj , Vj }j^z, acting on the sub- 



spaces Vj and Wj\ Aj : Wj 



Wj,Bj : Vj 



Wj,Tj 



Wj - 
as Aj 



Vj, where operators {Aj, Bj,Tj}j<=z are defined 
= Qjl Q,. I'., = (1,1 I ' r I. = P 3 TQ r The 



operator T admits a recursive definition via 



T- — 

1 3 ~ 



Aj+i B j+ i 



T 



where T, 



PjTPj and Tj works on Vj 



It should be noted that operator Aj describes interaction 
on the scale j independently from other scales, opera- 
tors Bj , Tj describe interaction between the scale j and all 
coarser scales, the operator Tj is an "averaged" version of 
Tj _ i . The operators Aj , Bj , Tj , Tj are represented by ma- 
trices a- 7 , fi> , 7-? , s- 7 



K (x, y)^j.k(x)^j, k ' (y)dxdy 



Tfe,fe' 



D k,k' 



K(x,y)ipj < k(x)(pj,k'(y)dxdy (4) 
K ( x > V^j.k (x)ipj,k> (y)dxdy 
K(x, y)<Pj,k( x )¥>3,k> (y)dxdy 



We may compute the non-standard representations of oper- 
ator d/dx in the wavelet bases by solving a small system 
of linear algebraical equations. So, we have for objects (4) 



CI; 



/ ip(2- j x - i)ip'{2- j - t)2- j dx 



PL = 



ill 



where 



2~ j J ip(2~ j x - i)(p'[2- j x - e)2- j dx 
J <p(2- j x - i)tp'(2- j x - £)2- j dx 

a e = 

h = 
it = 



ib(x — £)—ip(x)dx 
dx 

ijj(x — £)—ip(x)dx 
dx 

ip(x — £)—tp(x)dx 
dx 



then by using refinement equations we have in terms of fil- 
ters (h k ,g k ): 

L-l L-l 

Oij = 2^ g k 9k'r2i+k-k' , 

k=0 k'=0 
L-l L-l 

(3j = 2 ^ ^ 9khk'f2i+k-k' ; 

fc=0 fc'=0 
L-l L-l 

7i = 2^ hkg k >r 2 i+ k - k > , 
k=o fe'=o 

where re = j tp(x — £)4-(p(x)dx,£ G Therefore, the 
representation of d/dx is completely determined by the co- 
efficients or by representation of cZ/cfa; only on the sub- 
space Vq. The coefficients re,£ £ Z satisfy the following 
system of linear algebraical equations 



L/2 



|^>2fc-i( 



k=l 



-l\T2l-2k+l + r2e+2k-l) 



and Y,eJ r n = _1 > where 02fc-i = 2j]i = o hih i+2 k-i, 
k = 1, . . . , L/2 are the autocorrelation coefficients of the 
filter if. If a number of vanishing moments M > 2 then 
this linear system of equations has a unique solution with 
finite number of non-zero m, re ^ for — L + 2 < 



£ < L ~ 2, rg = — r_£. For the representation of op- 
erator d n /dx n we have the similar reduced linear system 
of equations. Then finally we have for action of operator 
Tj (Tj : Vj — ► Vj ) on sufficiently smooth function /: 

{Tjf)(x) = 2 ~ 3 ^2nfj,k-t I <Pj,k{x), 

kez \ e / 
where (fj } k(x) — 2-^ 2 ip(2- J x — k) is wavelet basis and 

f jt k-i = 2^' /2 / f(x)<p(2-ix -k + £)dx 



are wavelet coefficients. So, we have simple linear para- 
metrization of matrix representation of our differential op- 
erator in wavelet basis and of the action of this operator on 
arbitrary vector in our functional space. Then we may use 
such representation in all preceding sections. 
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